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ABSTRACT 


r/j  ftr*-' 

irn  and  main  beam  gain 


The  statistical  properties  of  the  power  pattern  and  main  beam  gain 
are  here  determined  for  an  array  of  randomly  located  submerged  acoustic 
sensors.  Tl.e  system  investigated  models  a  sparse  array  in  which  the  sen- 
cors  are  dropped  haphazardly  over  a  region  or  in  which  the  sensors,  how¬ 
ever  they  are  initially  placed,  become  spatially  diffused  by  a  process 
akin  to  a  two-dimensional  random  walk.  Signal  energy  is  assumed  to  arrive 
over  a  vertically  dispersive  channel  typical  of  the  long  range  deep  sea 
acoustic  channel  with  a  bigradient  sound  speed  profile.  Results  are  ob¬ 
tained  for  the  mean  value  of  the  power  pattern  and  the  mean  and  variance 
of  the  main  beam  power  gain  as  a  function  of  array  size.  It  is  shown  that 
in  typical  cases  the  Piean  array  gain  will  preserve  its  value  within  3  dB 
until  the  dispersion  parameter  of  the  array  measured  by. the  standard  dev¬ 
iation  of  element  location  is  around  35  wavelengths.  The  results  are  found 
to  be  consistent  with  results  obtained  by  others  for  the  coherence  distance 
in  a  multipath  acoustic  field. 
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STATISTICAL  PROPERTIES  OF  A  RANDOM  ARRAY  OF 
ACOUSTIC  SENSORS  IN  A  MULTIPATH  ENVIRONMENT 

INTRODUCTION 

The  statistical  properties  of  the  power  pattern  of  an  array  of  acous¬ 
tic  sensors  suspended  from  individual  freely  floating  buoys ,  receiving  from 
a  source  via  a  time-dispersive  medium  are  here  Investigated.  '  The  system 
geometry  is  Indicated  in  Figure  1.  The  array  elements,  numbering  N,  are 
assumed^  distributed  in  a  region  centered  on  the  origin  of  coordinates.  By 
independent  means  the  system  learns  the  position  (x^,  y^,  z^),  i  ■  1,2,  ...  N 
of  each  of  its  elements.  We  suppose  that  the  array  then  organizes  itself 
at  a  specified  frequency  forming  a  beam  aimed  in  some  selected  direction 
by  adding  suitably  phased  versions  of  the  element  outputs.  If,  for  Instance, 
the  element  outputs  were  of  unit  magnitude  at  the  specified  frequency,  if 
the  signal  were  arriving  at  angle  (6,$),  and  if  the  array  were  focused  to 
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X 


receive  a  signal  ray-  arriving  in  the  y-z  plane  at  a  colatitude  angle  6  the 

> 

complex  array  output  at  that  frequency  would  be, 


N  jk[*  sin6cos$  +  y  (sin9sin$ 

A-I*  ° 

n-1  r  " 


sin0  )  +  z  (coaf)  -  cos*  )] 

s  n  s 


(1) 


where  the  wavenumber  k  -  2tt/X*  X  being  the  wavelength.  In 
this  work  we  will  assume  that  all  elements  are  at  the  same  depth,  all  in 
the  x-y  plane,  so  that  all  zq  ■  0.  Though  surface  waves  will  cause  verti¬ 
cal  displacement  of  the  elements,  the.  wavelengths  of  interest  are  such  that 

« 

in  placid  seas  (sea  state  <  4)  the  displacement  is  less  than  0.1  wavelength. 
Furthermore,  the  system  envisioned  is  expected  to  utilize  a  drouge  with 
each  buoyed  element  so  that  vertical  motion  will  be  filtered.  Element  po¬ 
sitions  in  the  (x,y)  plane  will  be  assumed  independently  distributed  accord¬ 
ing  to  some  appropriate  two-dimentsional  probability  density  function  (pdf) 
f^(x,y),  as  will  be  discussed  later. 

The  dispersion  model  utilized  assumes  a  collection  of  M  planar  wave- 
fronts  impinging  on  the  array  all  originating  from  the  same  source  and 
all  arriving  with  the  same  azimuth  angle  but  with  different  colatitude  an¬ 
gles.  This  model  requires  that  the  phase  front  corresponding  to  a  given 
ray  arriving  at  the  array  center  be  adequately  approximated  by  a  planar 
surface  wherever  the  phase  front  contacts  the  array.  If  the  array  length 
In  the  azimuthal  direction  of  arrival  of  the  ray  is  d  and  the  colatitude 


angle  of  the  ray  is  6  the  distance 


across  the  wave  front  over  which  planarity 


should  hold  is  dcosO.  Each  wavefront  is  characterized  by  a  complex  amplitude 


B  e 
m 


m 


m  ■  1,2  ...  M , 


(2) 
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at  the  frequency,  f;the  phases  and  amplitudes  are  measured  at  the  origin 


of  coordinates,  the  former  relative  to  an  arbitrary 
be  treated  as  a  random  variable,  independent  of 
dom  variable  uniformly  distributed  in  2tt  and  $  for 

D 

sumed  independent.  .No  assumption  is.  made  about  the 
for  different  m.  The  colatitude  angle  of  arrival 


reference.  B  will 
m 

d>  will  be  assumed  a  ran- 
a 

different  m  will  be  as- 
dependence  among  the  B 

n 

\  \ 

of  a  wavefront  will  be 


denoted  0  ,  m  -  1,2,  ...  M.  The  0  will  be  viewed  as  nonrandom  constants. 

®  w  .... 

They  may  be  taken  to  be  equally  spaced  angular  samples.  The  modelj employed 
corresponds  to  one  used  by  Smith  [1]  to  calculate  spatial  coherence  in  a 
multipath  channel. 

With  the  wavefronts  arriving  at  azimuth  angle  $  and  the  array  focused 
to  receive  a  plane  wave  from  a  source  at  azimuth  angle  —•  ,  i.e.  from  a 
source  in  the  y-z  plane,  and  colatitude  angle  0fl  the  total  array  output  bi- 
comes 


M  N  jk[x  sin0  cos$  +  y  (sin6  sin$  -  sin0  )  +  ♦  ] 

. , .  ...  rr_^nm  n  m  s  m 

A(*,ea)  -  l  l  Bme 

This  is  the  conplex  array  pattern.  The  statistical  properties  of  the  cor- 

i  2  i  * 

responding  array  power  pattern  |A  ($,0g)|  will  be  investigated  below. 


The  terminology  used  here  differs  somewhat  from  that  given  in  Urick  [2], 

We  initially  calculate  an  array  pattern  |A^(<,0  )|  which  corresponds  to 
the  square  of  Uiick's  response  function  R^(0,$)  [2,pp.  49-50],  We  then  de¬ 
termine  a  normalized  mean  array  pattern  which  is  similar  to  Urick's  "beam 
pattern"  [2,  p.  50],  The  normalization  used  here  differs  however  from  that 
used  by  Urick  allowing  us  to  account  for  loss  of  coherence  across  the  array. 
Finally,  we  define  the  normalized  mean  array  pattern  evaluated  on  the  main 
beam  as  the  mean  power  gain. 

[1]  P.  W.  Smith,  Jr.,  "Spatial  Coherence  in  Multipath  or  Multimodal  Chan¬ 
nels, "  Jour.  Acoustic.  Society  of  America,  Vol,  60,  No.  2,  Augunt  1976, 
pp.  305-310. 

[2]  R.  J.  Urick,  Principles  of  Underwater  Sound,  McCraw  Hill  Book  Co., 
Second  Edition,  1975. 


STATISTICAL  PROPERTIES  OF  ARRAY  POWER  PATTERN 


The  mean  value  of  the  array  power  pattern  is  given  by 


M  M  M  H 


I  i  l<\\> 

Oj-1  m^-l  n^-1  1  2  > 

jktxQ  ein0ffl  cos*-xn  sin0n  cos4>+yn^(ain^ffl^8in(J)-aine8)-yn^(sih9B^8in(j)“Sin08) 


expectation  <Cexp(j*ffl  -  j*m  “  1  when  m^  "  m2  and  is  zero  other  ¬ 


wise  so  that  (A)  Is 


M  N  N  ,  jk(x  -x  )sin0  cos* 

<Ja2<M.)|>-  l  l  l  <>£><•  "i  2 

o-l  n,«l  nyl 


*1  2 


Jk(yn  -yn  )  (sin0msin<ji  -  siut)#)„> 
u  e  12 

(5) 

Using  the  assumption  that  the  element  positions  are  independent  random  vec¬ 


tors  (5)  is  written 


M  N  N  jk(x  -x  )sin0  cos^ 

<_|A2(*.8s)  l>-  l  <?l>  t"  +  I  !.<•  T2  .. 

n,-l  n«"i 


1  A  ‘2 

"l’S 


i <»:>!*! 


m  > 

0-1  ln-i 


jk(y  **y  )  (sinO^sin*  -  sln08)^> 

•  e  nl  n2 

.2 

.  N  Jktx  sin0  cos*  +  y  (slnO  sin*  -  sin0g)l> 
r  -  n  n  n  m  1 

+  !<e  1 


'  N  Jk[x  sin8  cos*  +  y  (sin0  sin*  -  sinO  )]  , 

-  I  n  m  s  ->{2  (6) 


-5- 


The  expectations  inside  the  brackets  are  two  dimensional  characteristic (func¬ 
tions  of  tht?  random  vectors  (xn,yft).  Assuming  all  elements  to  have  identi¬ 
cally  distributed  location  vectors  and  denoting 


jk[x  sln0  cos$  +  y  (sin(J  sin$  -  sin0  I 

<«  "  I  <7> 


then 


M 


CIA  <+*0a)l>-  I  <b;>[n +(NZ-N)|« xy<en.*.e8)h 

m«l  1 


(8) 


At  this  point  we  specialize  the  distribution  of  the  location  vectors. 
We  assume  the  effect  of  the  forces  tending  to  scatter  the  array  elements 
to  be  modeled  by  a  two  dimensional  random  walk  with  independent  increments 
along  the  coordinate  axes.  (xn,yn>  aff®r  a  time,  be  distributed  ac¬ 


cording  to  a  two  dimensional  random  variable  approaching  a  normal  withj  vari- 

2  2 

ance  along  x  and  y  given  by  a  and  c  respectively.  In  this  case 

x  y 


*  (0  )  * 

xy v  m*  '  s'  , 


-[02k2sin2e  cos2<J>  +  a2k2(sin0  sin<fr  -  sin0s)2]/2 
e  x  m  y  a  , 


(9) 


The  model  chosen  accounts  for  element  diffusion  but  Ignores  drift  compon¬ 
ents  which  are  sure  to  be  present.  The  assumption  is  Implied  that  trans¬ 
lation  of  the  entire  .<rray  will  not  seriously  affect  array  response  when 

attempting  to  focus  on  a  distant  target  if  the  translation  is  small.  Wc 
1"  2 

point  out  that  a  and  o  are  functions  of  time;  as  the  array  ages  these 
x  y 

parameters,  which  are  a  measure  of  the  size  of  the  array,  will  grow. 
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If  we  were  to  suppose  that  initially  the  elements  are  close  to  one 
another  so  that 


♦  (e  ,$,e  )  -  l 

xyv  m*  *  s 


<|A2a,0  )1>-  N2  l  <B2>  -  N2B2  (11) 

*  o-l  “ 

B2  -  J  ^  B2  ^  Is  the  mean  square  value  of  the  total  signal  arriving  at  the 
array  center  and  (11)  represents  the  power  delivered  by  the  array  when  it 
is  sufficiently  small  not  to  be  defocused  by  the  multipath.'  It  is  useful 

to  normalize  |a2(4>,6  )|  by  N2B2.  We  denote  this  random  variable/ 

8 


r  -  — 


I  a2 (♦ *  e  )  | 


the  normalized  array  power  gain.  Its  mean  value  is 


1  .  n  1,  ?  <\>  “V  8ln V08  *  -02k2(sinensin*-slnesy 

<  ^  U  b2  6  (13) 


*An  alternative  normalization  is  given  by 

|a2om8)! 

ri  "  .  ~H  j*  o 

-2I  l ».•  "I2 

m-1 

tbl.  random  variable 

to  small  size.  T  as  given  by  (12;  PP  total  arriving  signal 

is  unity,  in  r,  an!  the 

power  at  the  array  center  has  been  1 {c™v*  *  £f  t  of  prray  defocusing. 

quantity  is  more  nearly  representative  of  .  w  h ave  there- 

?,  is  however  e  more  difficult  ,usntity  with  which  to  work  We  have  tn 
fore  settled  on  t  in  which  the  normalisation  is  done  with  a  constant. 
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As  a  final  step  we  define 


0(6  )A6 
n 


(U) 


that  Is,  we  suppose  the  continuum  of  possible  signal  arrival  angles  to  be 

quantized  Into  increments  48  and  that  the  fractional  power  obtained  from 

« 

the  mth  Increment  is  0 (8^)46.  By  allowing  the  Increments  to  become  small 
(13)  will  be  approximated  by  ar.  Integral  as  follows: 

,  ,  -oL\^8in^Qcos^<|>  -0^1t^(sin9sin<j>-sin9  )* 

<T>  -  ~  +  (!-|)  /  0(6)e  X  y  8  d9 

i  (15) 


<r>  as  given  by  (15)  is  shown  evaluated  in  Figure  2  for  the  case 

# 

2  2  2  I 

o*  m  or  -  0  with  ok/2fr  -  o/X,  a  family  parameter,  given  by  5,10,20,  and  40  wave 

*  J  I 

lengths.  The  power  density  0(0)  is  assumed  uniformly  distributed  over  ±10* 
relative  to  the  horizontal.  The  array  is  assumed  focused  for  a  source  in 
the  plane  of  the  array, ^that  is  9g  ■  90°,  and  at  an  azimuth  angle  of  90°. 

Of  particular  interest  is  the  magnitude  of  the  normalized  mean  array 
pattern  evaluated  on  the  main  beam  as  a  function  of  array  size.  We  refer 
to  this  quantity  as  the  power  gain  I*  .  Its  mean  value,  is  obtained 

from  (13)  or  (15)  evaluated  at  $  "  90* *  Thus  using  the  discrete  ray  model 


(16) 


FIGURE  2 


.  PATTERN  OF  NORMALIZED  MEAN 


ARRAY 

POWER  GAIN. 

(a) 

o/A  * 

5 

(b) 

o/A  “ 

10 

(c) 

o/A  " 

20 

<d) 

o/A  - 

40 

I 

Or,  using  the  continuous  approximation  for  the  Integral, 

,  ,  -oV  (sine-sine.)2 

<'rd>"  N  +  <1N>  /  0(0>a  d8  (17> 

la  (16)  and  (17)  the  subscript  on  a2  has  been  dropped.  (17)  has  been 
evaluated  numerically  as  a  function  of  the  normalized  size  variable 

0k/2v  «  o/X,  with  the  angular  distribution  of  energy  arriving,  6(6),  aa  a 
parameter  and  6*  set  to  zero.\  8(8)  was  assumed  uniform  over  angles 
is*,  ilO*,  and  +20%  relative  to  the  plane  of  the  array.  Results  are  shown 

in  Figure  3,  For  distant  sources  the  arrival  angles  are  apt  to  be  within 

*  ^ 

NORMALIZED  MEAN  POWER  GAIN 


FIGURE  3.  j  MEAN  POWER  GAIN  AS  A  FUNCTION  OF  NORMALIZED 
}  ARRAY  SIZE  VARIABLE,  a/X. 

(a)  8(6)  -  1/10,  85*  <6<  95* 

(b)  8(6)  -  1/20,  80*  <B<  100’ 

(c)  8(6)  -  1/A0,  70*  <e<  no* 


•M*  ■ 


ilO*.  Note  that  delivered  power  is  reduced  to  about  1/2  when  o/X  *  35  wave¬ 
lengths  J  For  sources  nearby,  bottom  and  top  reflections  may  result  in 
energy  arriving  at  steeper  angles  and  the  +20%  distribution  may  be  viewed 
as  a  model  suggesting  the  effect  in  such  a  case.  Here  delivered  power  is 
reduced  to  about  1/2  where  o/X  *  8  wavelengths. 

Setting  the  angle  0g  to  zero  means  focusing  the  array  for.  signal  ar¬ 
rivals  in  the  plane  of  the  array.  This  is  not  optimum  for  signals  arriv¬ 
ing  over  a  dispersion  of  latitude  angles.  To  show  this  we  hare  plotted 
this  mean  power  gain  as  a  function  of  the  colatitude  angle  0  for  the  case 
of  8 (0)  uniform  in  +10*  around  the  horizontal  and  o/X  -  33.3  wavelengths. 

The  result  is  shown  in  Figure  4.  The  maximum  is  seen  to  occur  with  the 

beam  formed  for  0  -  86*  (by  symmetry  it  will  also  be  maximum  for  0  ■  94#). 

*  8 


80  82  84  88  88  90 

COLATITUDE  ANGLE,  &s 


FIGURE  4.  MEAN  POWER  GAIN  AS  A  FUNCTION  OF 
VERTICAL  AIMING  ANGLE  0g. 

B(©>  -  1/20,  80*  <9<  100° 
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The  variance  of>  the  power  gain  will  be  useful  as  an  indicator  of  the 

2  v.  2 

gain  variability.  We  thus  determine  Var  -  (<To>  -  <1  Tj-  )  where 
T  «  | A2  (■?■,  6  )  |/N2B2.  Starting  (3)  evaluated  at  $  “  w/2  we  have 

O  &  8 


/(2*0s>  2  Ao  "  l  I  lBn' 
0*1  n*l| 


M  N  j[kyn(sin0m  -  sinS^  +  4>01 


The  fourth  moment  of  the  magnitude  is  given  by 


C|A4|>>-  I  l  <B  B  B  B 

^ '  o  4  4  m,  nu  m,  m, 

n^*n2,n3»  12  3  4 


o4  -  1  n4  «  1 


_ 1  (k.  y  -k  y  +k  y  -k  y 

.  <C  nl  nl  "2  “2  m3  n3  “4  n4 


.<eJ<*«1'S+S'*"4>>  <19> 


We  have  used  the  abbreviated  notation  k  ■  k(sin0  -sin0  ).  With  the  6 

m^  m^  s  ifl^ 

independent  uniformly  distributed  random  variables  in  (0,2tt)1(19)  reduces  tc 


«  -  M  j[k  (y  -y  )+k  Cy  -y  )] _ 

<'<>  -  ,i  i  1  1  2  3  3  4  ^ 


^1 ^n3 


l .  <b„2>  <*J> 

VV1  1  2 

n»2 


J  [k  (y  -y  )  +  k  (y  -y  )  3 

fl«2 


u  3(k  (yrt  -y  -y„ 

M  ,  m  'n,  n.,  n.  ~j> 

*  l  <.»!><•  1  34 

m-1 


waafr 
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Suomarion  over  the  n^,  1  *»  1,2,3,4,1s  now  carried  out  using  the  assump¬ 
tion  that  the  positions  of  the  different  elements,  the  y^  k-1,2,  ...H, 
are  independent.  The  summation  is  straightforward  although  laborious. 
Carrying  out  the  nteps  we  can  then  write 


Var  T  - 
o 


<ia*i>-  <ia;i>2 

n4b4 


M  m  <■£><£>  r  r-A*  -A*n 

-!  !  — 1  -4-'-2  ■-\  W(h-dL  x+. 

tiij“l  c\2**l  B  N  ^ 

r-o2(k  +k  )2  -02(k  -k 

®i  °1  m2 

+  2N(N-1)  |e  1  1  +  e  J 

-V(k2 +k2  )r-T-<k  +k  >2  ~T^k  -k  )21 

2 ->  •  m1  -  m9  2  m1  m-  2  m.  nu 

+  4N(N-l)(N-2)o  +  e  J 


[2N(N-1) (N-2)(N-3)  -  NZ(N-1)  ;  e 


-o2 (k2  +k2  )\ 

mi  m2  r 


+  £  JZ-3 - - - E -  .  -L  ^N-l)  +  4N(N-1)  e 

_ i 


-o2k2 
2  0  m 


+  N(N-l) (N-2) (N-3)e 


-2c2k2 

m 


+  2N(N-1) fN-2)  e 


-3o2k2 

m 


+  N(N-l)  e 


-4o2k2] 


The  ratio  ’’ >r  T  i<^X  >  ,  where  the  numerator  is  given  by  (21)  and 
o  o 

the  denominator  <.y  the  square  of  the  mean  power  gain  given  in  (16),  is 
a  useful  measure  of  relative  variance.  The  result  is  cumbersome. 


K  -  'Vi*/. 
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however,  and  expressions  applicable  to  limiting  cases  are  instructive. 

t 

Tv?o  such  cases  are  here  evaluated  assuming  the  ray  amplitudes,  B  ,  are 

k  m 

equal  for  all  m  and  constant  so  that<C.B  ,/B  *  (1/M)2,  k  even. 

m  *’ 

The  first  case  treated  assumes  the  elements  very  widely  dispersed 
2  2 

so  that  Ok  »  1  for  all  rays  except  ones  for  which  sinQ  *  sii;9  . 

m  ms 

Since  k  *■  k  (sin9  -  sin6  )  a  ray  along  the  aiming  angle  will  result  in 
ni  m  8 

V.  ■  0.  Assuming  one  ray  is  along  the  aiming  angle  set  at  0g  •  90° 
we  get,  using  (16)  and  (21), 


Var  T 
_ c 

<ro>; 


NM2  +  M(2N2-2N-1)  ~  (2N2-H-1) 
2 


N(M+N-1) 


(22) 


When  M  "  1,  there  is  only  1  ray,  and  that  one  along  the  aiming  angle. 
There  is  no  multipath  and  the  array  will  be  correctly  focused.  The 
ratio  above  is  then  zero.  When  M  gets  large  without  bound  while  N  re- 
mains  finite  the  ratio  approaches  unity.  This  result  can  be  anticipated. 
For  M  large  the  elements,  being  widely  dispersed,  see  a  sinusoid  with 
Rayleigh  magnitude  and  random  phase.  The  ratio  in  (22)  in  that  case  is 
that  of  the  variance  and  squared  mean  of  an  exponential  random  variable 
for  which  this  ratio  is  unity.  If  M  is  held  finite  while  N  is  allowed 
to  increase  without  bound  the  ratio  tends  to  zero.  This  result  arises 
because  the  element  outputs  caused  by  the  one  ray  along  the  aimiug 
angle  are  coherently  combined  by  the  array.  The  random  component  con¬ 
tributed  by  the  rays  off-axis  add  up  non-cbherently  at  the  array  output. 
The  latter  are  the  variance  producing  components.  But  as  N  Increases 


This  result  assumes  also  that  sinkra  ^sinkm  for  all  Should  there 

12  x  z 

be  rays  arriving  symmetrically  relative  to  the  horizontal  there  will  be 
some  sink,,,  “sinkj^  for  m./m?.  In  such  a  case  some  additional  terms  will 
be  required  from  ^21)  and  (22)  may  exceed  unity. 
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The  ratio  of  the  non-coherent  components  tends  to  zero. 

We  point  out  that,  if  no  ray  comes  in  at  an  angle  sufficiently  close 

2 

to  0g  to  make  a  significant  coherent  contribution  then  Var  F0/<Lr'o^> 
approaches  (l-l/MN).  Now  the  ratio  approaches  unity  with  increasing 
M  or  N  as  one  would  expect. 


The  second  case  treated  is  one  for  which  the  element  locations  are 

reasonably  compact,  say  within  a/\  ■  10,  and  N  is  large.  In  this  case, 

2  2 

if  ray  arrivals  are  within  i  10*  of  the  vertical  focusing  angle  6  ,  o  k 
-to***  •  ■ 

is  small  and  e  is  close  to  unity  (it  is  0.748  for  0  *  10°).  The 

predominant  terms  in  (21)  are  then  those  of  highest  degree  in  N.  Extract- 


i 


In  (23)  we  have  set  such  factors  as  (N-l) (N-2) (N-3)/N^  to  unity  in  con¬ 
formity  with  specification  that  N  be  large.  Again  now  with  a  constant 
and  equal  for  all  m 
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2 

is  obtained  from  (16)  using  the  assumptions  pertinent  to  this 

case.  The  ratio  inside  the  bracketed  facto,  in  (24)  is  1/M  for  o  «  0. 

For  a  “  0  (24)  is  identically  the  variance  of  the  squared  amplitude  of 

the  sum  of  M  equal  amplitude  sinusoids  with  independent  phases,  all 

uniformly  distributed  in  2tt;  this  result  for  a  ■  0  simply  reflects 

the  variability  of  the  incoming  total  signal  magnitude  at  the  array 

center.  It  is  interesting  to  observe  that  if  B  were  Rayleigh  dis- 

m 

tributed,  Var  would  be  unity  for  all  M  and  for  all  o/X 

for  which  (23)  is  valid. 

To  convey  some  idea  of  how  the  variance  changes  with  the  array  size 

parameter  a/X,  and  to  provide  some  results  without  approximations  the 

2 

ratio  Var  ro/<^* To was  calculated  for  some  representative  cases 

using  (16)  and  (21)  as  they  stand.  Table  1  shows  these  results  for  the 

case  of  M  rays,  M  ■  3,  5,  11,  and  21,  the  rays  arriving  at  equally 

spaced  angles  in  an  interval  of  £  10°  relative  to  the  horizontal.  The 

beam  pointing  angle  was  set  along  the  horizcntal  (0fi  ■  90°)  and  in  the 

direction  of  the  source.  The  number  of  sensors  N  was  taken  to  be  31. 

It  is  worth  noting  that  for  M  *  3,  N  ■  31  in  the  limiting  expression 

(22)  (o/X-*»),  Var  V  ^ ^  •  0.1175  only  slightly  below  the  calculat* 

o  o 

ed  value  in  Table  1  for  M  ■  3,  and  u/X  ■  60. 

There  may  be  another  phenomenon  responsible  for  the  trend  of  the 
calculated  values  in  Table  1.  With  M  >  1  the  signal  power  at  the  array 
center  is  a  random  variable  as  a  consequence  of  the  random  phases  of  the 
incoming  rays.  The  ratio  of  Var  rQ/<^ro^>  is,  as  pointed  out  below 
(24),  given  by  (1-1/M).  This  is  exactly  the  value  in  Table  1  for 


>*1Q 


o/A  ■  0.  As  at  A  increases  the  ratio  Var  rQ/<^  To  >  decreases,  pertly 
for  the  reason  discussed  under  (24) ,  and  perhaps  also  because  the  signa 
amplitude  which  varies  with  position  is  averaged  by  the  array. 


M 


3 

5 

11 

21 

0 

0.67 

0.80 

0.91 

0.95 

5 

0.66 

0.80 

0.91 

0.95 

10 

0.62 

0.78 

0.91 

0.95 

15 

0.42 

0.73 

0.90 

0.95 

20 

0.20 

0.70 

0.88 

0.95 

40 

0.12 

0.64 

0.84 

0.92 

60 

0.12 

0.47 

0.81 

0.91 

TABLE  1.  RATIO  OF  VARIANCE  TO  SQUARED  MEAN  OF  ARRAY 
OUTPUT  AS  A  FUNCTION  OF  ARRAY  SUE,  o/A,  AND  NUMBER, 
M,  OF  EQUAL  AMPLITUDE  RAYS.  RAYS  ARE  ASSUMED  EQUALLY 
SPACED  IN  AN  INTERVAL  ±  10  FROM  HORIZONTAL. 
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CONCLUSION 

Mean  power  pattern,  mean  main  beam  gain  and  main  beam  gain  variance 
have  been  determined  for  an  acoustic  array  of  widely  scattered  submerged 
elements  The  elements  are  assumed  to  be  organized  to  accept  a  plane 
wave  but  the  array  sees  a  multipath  field  typical  of  the  acoustic  field 
at  a  great  distance  from  a  source  in  the  deep  ocean.  The  results  obtained 
show  the  diminution  of  array  effectiveness  as  the  array  size  grows;  for 
a  typical  case  where  the  range  of  latitude  angles  of  the  arriving  signal 
is  10*  above  and  below  the  horizontal,  the  mean  power  gain  fal.1  3  dB  when 
the  element  spread  as  measured  by  the  element  position  standard  devia¬ 
tion  is  35  wavelengths. 

In  earlier  work  Smith  [1]  calculated  the  normalized  coherence  magni¬ 
tude  as  a  function  of  horizontal  separation  in  a  long-range  transmission 
channel  with  bigradient  sound  speed  profile.  In  particular  he  explicitly 
obtains  the  coherence  distance  for  50%  coherence  for  a  source  on  axis  and 
a  receiver  close  to  the  channel  edge.  It  turns  out  to  be  46  wavelengths 
if  the  difference  in  sound  speed  between  the  axis  and  the  receiver  loca¬ 
tion  is  20m/sec.  The  range  of  vertical  angles  of  arrival  is  ±  9.36*  rela¬ 
tive  to  the  horizontal  at  the  receiver  and  the  energy  density,  obtained 
from  earlier  results,  is  assumed  by  him  to  be  uniform  over  the  range  of 
arrival  angles.  This  situation  is  roughly  the  same  as  that  used  to  ar¬ 
rive  at  the  35  wavelength  figure  mentioned  above.  The  two  results  are 
interestingly  similar.  This  ought  to  come  as  no  surprise;  since  the  co¬ 
herence  distance  corresponds  to  the  array  dimension  useful  for  coherent 
combination  of  the  spatially  sampled  field.  The  crucial  factor  in  de¬ 
termining  the  coherence  distance  is  the  range  of  angles  of  arrival.  To 
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see  this  Intuitively,  imagine  two  multipath  rays,  one  horizontal  and  the 
other  smaller  in  amplitude  and  at  9*  to  the  horizontal.  The  phase  of 
the  resultant  vs  horizontal  distance,  x,  is  expressible  as  [3] 

*Kx)  -  kx  +  $  (x) 

where  k  is  the  wave  number.  $(x),  is  periodic  with  period  2ir/k(l‘cos0)  ■ 
X/(l-cos6).  With  0  ■  10*  the  period  is  about  66X.  This  is  not  too  far 
different  from  the  2o  range  (-  70  wavelengths)  for  an  array  power  gain  of 
0.5,  or  from  the  46  wavelengths  obtained  by  Smith  for  coherence  distance 
for  0.5  normalized  coherence. 


[31  F.  Haber,  ’’Phase  Variations  with  Position  in  an,  Underwater  Multi' 
path  Environment  and  its  Effect  on  Array  Pattern,"  Valley  Forge 
Research  Center  Quarterly  Progress  Report  No.  24,  University  of 
Pennsylvania,  the  Moore  School  of  Electrical  Engineering,  Phila¬ 
delphia,  PA  19104,  pp.  20-30. 


-19- 


ACKNOWLEDGMENTS 

I  am  Indebted  to  Mr.  William  J.  Graham  for  hla  review  of  part  of 
this  work  and  for  his  suggestions  for  improvement,  to  Mr.  DeYuan  Ho  for 
hie  help  with  .he  numerical  evaluations,  and  to  the  U.  S.  Office  of 
Naval  Research  for  its  support  on  Contract  N00014-77-C-0252.  X  am  also 

indebted  to  the  personnel  of  the  Office  of  Naval  Research  and  of  the 

« 

Naval  Ocean  Systems  Center,  San  Diego,  California,  for  providing  a 
system  description  which  was  used  as  a  guide  to  practical  values  in 
this  paper. 


-M<«  I- 


•|H»  III  ~t  "I  i 


xmaa 


na 


I 


/~ 


11  '"’(U)  RESEARCH  IN  DISTRIBUTED  ARRAYS 

1  -  AGENCY  ACCESSION  NO-  DN77539© 

5  -  DATE  OF  SUNNARY  ■  12  NOV  30 

—  39  -  PROCESSING  DATE  (RANGE)-  30  NOV  SO 

6  -  SECURITY  OF  WORK-  UNCLASSIFIED 

—  12  -  S  *  T  AREAS- 

oooaoo  ACOUSTICS 
000100  ACOUSTIC  DETECTION 

—  aiE  -  N1LITARY/CIVILIAN  APPLICATIONS'  MILITARY 

— 10A1  -  PR 1 NARY  PRDGRAN  ELENENT'  S2711N 

— 10A2  -  PRINARY  PROJECT  NUMBER •  F11121 

--10A2A  -  HRINARY  PROJECT  AGENCY  AND  PRDGRAN •  Filial 

— 10A3  -  PRINARY  TASK  AREA'  RF11 131810 

--10A4  -  UDRK  UNIT  NUMBER •  NR-139-102 

--17A1  -  CONTRACT, “GRANT  EFFECTIVE  ..^E*-  MAR  7? 

--17A3  -  CONTRACT, “GRANT  EXPIRATION  DATE'  FEB  30 

—  17B  -  CDNTRAC T /GRANT  NUNBER'  A0©014-??-C-©25a 

—  17C  -  CONTRACT  TYPE'  COST  TV) E 

--17D3  -  C ONTR AC T /GRANT  ANOUNT •  S  25,000 

--  17E  -  KIND  OF  AWARD •  EXT 

--  17F  -  CONTRACT  “GRANT  CUNULAT1VE  DOLLAR  TOTAL*  >  130,000 

—  19A  -  DOD  ORGANIZATION-  OFFICE  OF  NAVAL  RESEARCH  322 

1  OF  1 

—  11  -  TITLE*  (U)  RESEARCH  IN  DISTRIBUTED  ARRAYS 

—  1  -  AGENCY  ACCESSION  NO-  DN775a9© 

—  a  -  DATE  OF  SUNNARY'  12  NOV  SO 

—  39  -  PROCESSING  DATE  (RANGE)-  3©  NOV  30 

6  -  SECURITY  OF  UORK'  UNCLASSIFIED 

—  13  -  S  T  AREAS- 

oooaoo  ACOUSTICS 
000100  ACOUSTIC  DETECTION 

—  aiE  -  NILITARY /CIVILIAN  APPLICATIONS'  DILI'  iRV 

— 10A1  -  PRINARY  PROGRAM  ELEMENT •  82711N 

— 10A3  -  PRIMARY  PROJECT  NUMBER-  Filial 

--lOAaA  -  PRINARY  PROJECT  AGENCY  AND  PROGRAM'  Filial 

— 10A3  -  PRIMARY  TASK  AREA'  RF1 1131810 

--10A4  -  UORK  UNIT  NUMBER'  NR -129- 103 

— 17A1  -  CONTRACT/GRANT  EFFECTIVE  DATE-  MAR  77 

— 17A3  -  CONTRACT/GRANT  EXPIRATION  DATE-  FEB  SO 

—  17B  -  CONTRACT/GRANT  NUMBER'  N00014-77-C-0253 

—  17C  -  CONTRACT  TYPE-  COST  TYPE 

— 17D3  -  CONTRACT/GRANT  AMOUNT •  »  25,000 

—  17E  -  KIND  OF  AWARD-  EXT 

—  17F  -  CONTRACT/GRANT  CUMULATIVE  DOLLAR  TOTAL*  I  130,000 

—  19A  -  DOD  ORGANIZATION-  OFFICE  OF  NAVAL  RESEARCH  233 


—  19B  -  DOD  ORG.  ADDRESS'  ARLINGTON.  VA  23317 

—  19C  -  RESPONSIBLE  INDIVIDUAL'  BOYER.  G  L  233 

—  19D  -  RESPONSIBLE  INDIVIDUAL  PHONE-  303-898-4305 

—  19U  -  DOD  ORGANIZATION  LOCATION  CODE*  5110 

—  19S  -  DOD  ORGANIZATION  SORT  CODE'  35333 

—  19T  -  DOD  ORGANIZATION  CODE •  265350 


20A  -  PERFORMING  ORGANIZATION'  L'NIV  OF  PENNSYLVANIA.  MOORE  SCHOOL  OF 
ELECTRICAL  ENGINEERING 

SOB  -  PERFORMING  ORG  ADDRESS •  PHILADELPHIA,  PA  19174 
ROC  -  PRINCIPAL  INVESTIGATOR'  HABER,  F 

ROD  -  PRINCIPAL  INVESTIGATOR  PHONE*  215-243-3104 
EOF  -  ASSOCIATE  INVESTIGATOR  (1ST)'  HABER,  F 
ROG  -  ASSOCIATE  INVESTIGATOR  ( END >  •  LIN,  T 
SOU  -  PERFORMING  ORGANIZATION  LOCATION  CODE •  4201 
RON  -  PERF.  ORGANIZATION  TYPE  CODE-  1 
20S  -  PERFORMING  ORG .  SORT  CODE*  31456 
EOT  -  PERFORMING  ORGANIZATION  CODE'  237000 

E2  ~  KEYWORDS'  (U>  RANDOM  ARRAYS  /(U>  DISTRIBUTED  ARRAYS/ (U > 

SONOBUDY  /(U)  PASSIVE  SONAR  / 

37  -  DESCRIPTORS'  Ui>  DISTRIBUTION /<U>  GAIN  /(U>  PASSIVE  SONAR 
tU)  SONOBUOYS  /(U>  UNDERWATER  SOUND  ;<U)  ACOUSTIC  COMPATIBILITY 
(U>  ^UNDERWATER  ACOUSTICS  /(U)  I ACOUSTIC  ARRAYS  /(U>  ^ACOUSTIC 
DETECTION/ tU>  ACOUSTIC  TRACKING  /<U>  EFFECTIVENESS  /(U> 

ELECTROMAGNETIC  COMPATIBILITY  / 

S3  -  TECHNICAL  OBJECTIVE'  <U)  TO  DEVELOP  TECHNIQUES  IN  SIGNAL  ANALYAIA 
TO  MAXIMIZE  ACHIEVABLE  ARRAY  GAIN  OF  A  RANDOM  SONOBUOY  ARRAY. 

24  -  APPROACH'  CU>  INVESTIGATE  PHASE  DECORRELATION  EFFECT  ON  ARRAY  GAIN 
AND  WAYS  TO  OVERCOME  THESE  EFFECTS. 

25  -  PROGRESS'  <U)  EARLIER  WORK  FOCUSED  ON  METHODS  OF  LOCALIZING  ELEMENT 
OF  A  RANDOM  ARRAY.  THEORY  HAS  BEEN  DEVELOPED.  VALLEY  FORGE  RESEARCH 
CENTER  QUARTERLY  PROGRESS  REPORT,  FEB.  1979  (U). 

13  -  WORK  UNIT  START  DATE-  MAR  77 

14  -  ESTIMATED  COMPLETION  DATE'  CD 


